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We theoretically investigate the role of different phases of coupling constants in the dynamics of atoms and 
two cavity modes, observing deterministic generation of prototype or hybrid Bell, W, GHZ, and cluster states. 
Commonly induced dipole-dipole interactions (far-off resonance) are inhibited between particular pairs of qubits 
under suitable choice of those phases. We evaluate the generation fidelities when imperfections such as dissipa- 
tive environments and time precision errors are considered. We show violation of local realism for the generated 
cluster state under such imperfections, even when approaching the weak coupling regime. 
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I. INTRODUCTION 

Near to complete half a century since its first description, the 
Jaynes-Cummings Model (JCM) [1] still remains at the core 
of outstanding research. Most of the advances in new tech- 
nologies and devices for quantum computing implementation 
based in solid state or atomic and optical systems made use 
of this model. Many tasks in quantum computing demand co- 
herent superpositions and/or entanglement. These two intrinsi- 
cally quantum ingredients are thus fundamental to disclose the 
quantum content of information. There has been a great effort 
to successfully and faithfully implement quantum computing 
building blocks in large scales. The simplest block is an in- 
formation carrier capable of being in general superposition of 
two orthogonal states, also known as qubit [2]. Another fun- 
damental block is the ability to perform general quantum state 
manipulations and measurements on these qubits [2]. 

The increasing degree of control that has been achieved in 
the level of single quantum systems has confirmed the predic- 
tions of fundamental models like the JCM and some of its gen- 
eralizations [3]. The inclusion of a second radiation mode (bi- 
modal cavities) is expected to bring out different quantum ef- 
fects [4—7]. Such effects may concern quantum statistical prop- 
erties of the fields, and also the dynamics with applications in 
the controlled generation of entanglement. The interaction of 
two-level atoms with two modes sustained by two weakly cou- 
pled cavities separated by a thin dielectric barrier is considered 
in [8], where the authors propose a generation scheme of atom- 
atom maximally entangled states. From the experimental side, 
the preparation of two modes of a superconducting cavity in a 
maximally entangled state has already been achieved [9]. Fol- 
lowing this experimental achievement, generation schemes of 
important multipartite entangled states such as GHZ, W and 
cluster states [10, 11] have been suggested. In general, such 
schemes make use of a sudden change of the atomic transi- 
tion frequency in order to induce resonant interaction with one 
mode or the other This might introduce imperfections [12] 
that should be considered for the best agreement with the ex- 
periment. 

In this paper, we study the dynamics of two non-degenerate 
bosonic modes interacting with N identical two-level systems 
(qubits), aiming at the generation of important multipartite en- 
tangled states between the modes, qubits or modes and qubits 



(hybrid). An important feature of our proposal is that it does 
not make use of sudden frequency changes. By fixing the qubit 
frequency separation to be midway the natural frequencies of 
the modes [8], and by properly choosing the phases of the cou- 
pling constants, our schemes avoid the undesired effects in- 
troduced by suddenly switching the frequencies. Besides the 
generation of Bell, W and GHZ states, we also describe a 
one-step generation scheme of the linear cluster state among 
four qubits [13]. When comparing, for example, with cavity 
QED schemes, it does not make use of either rotation of qubits 
(Ramsey zones) or classical external fields, usual requirements 
in most schemes for entanglement generation in these setups 
[14, 15]. It is remarkable that the effective model obtained 
here (XY model) allows the generation of such a state, com- 
plementing the usual approach which employs the Ising model 
[13]. It is important to emphasize that all protocols for the 
generation of multipartite entangled states presented here are 
deterministic, i.e. no measurements are performed at all. 

We organize the paper as follows. In Section II, we describe 
the model and its dynamics in the case of approximate reso- 
nance of the qubits with the bosonic modes considering: (i) 
not simultaneous interaction of qubits and modes, and (ii) si- 
multaneous interaction of qubits and modes. We then discuss 
the generation of Bell and W states. In Section III, we derive 
the effective Hamiltonian of the qubits simultaneously inter- 
acting with the modes in the case of far off resonance, and we 
propose the generation of Bell, W, GHZ and cluster states in 
this configuration. In Section IV, we consider common exper- 
imental imperfections such as amplitude damping and finite 
time-control precision, analyzing the fidelity in the generation 
schemes and the violation of local realism for the cluster state. 
In Section V, we present physical setups where our ideas may 
be implemented. We summarize and conclude in Section VI. 
The ideas in this paper will be put in the general language of 
two-level systems (qubits) and bosonic modes. Consequently, 
the results presented here may be of interest for application in 
any particular quantum system described by JCM interactions. 
This includes circuit-QED systems [16] and trapped ions [17], 
just to name a few. 



II. GENERATION OF W STATES WITH BIMODAL JCM 
INTERACTIONS 
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In cavity-QED, phase differences between coupling con- 
stants of atoms and field modes are in general considered 
unimportant in many situations. However, special choices of 
these phases have been shown to avoid the destruction of mi- 



croscopic correlations between the modes, important for the 
appearance of decoherence-free subspaces and robust states 
[18, 19]. Besides, it can not be underrated that these phases 
are of key importance for a plethora of quantum interference 
effects, responsible for many quantum phenomena [20]. In the 
context of trapped ions, adjustments of intensity and phase of 
the couplings constants have also been proved to be fundamen- 
tal to the generation of cluster states when the ions are driven 
by red- and blue-sideband laser fields [21]. 

In order to explore such possibilities, we study a general 
system composed of N two-level systems or qubits interacting 
with two bosonic modes. In the rotating wave approximation, 
the Hamiltonian for the two modes and qubit k {k = 1 . . . N), 
in the Schrodinger picture, reads (fi = 1) 



with coefficients 



,(S) 



n'^\b^a^+bat), 



(1) 



where lua (^b) is the frequency of mode A (B), a and b (a^ and 
b^) are the photon annihilation (creation) operators of modes 
A and B, lo is the qubit frequency and a], (i — x, y, z), with 
CT^ — (1/2) (cr^ ± icr^), refer to the Pauli operators of qubit k. 
Here we consider the qubit coupling with mode A to be real, 
positive, and constant i.e. Qj^ ' = fl for all k, while the qubit 

coupling with mode B may also be negative, i.e. Qj^ ' ~ fisfe, 
where Sk = ±1. 

We move to a frame rotating according to e^**^*, with 

G = LJa}a + ujb'^b+-(jl, (2) 

so that the Hamiltonian in this frame is given by 

Hk = Aa^a - A6^5 + ri(aV^ + aa+) + Osfe(6V^: + 6cr+), 

(3) 
where we have set the qubit frequency midway the modes 
frequencies separation, i.e. the detuning is given by A = 

WA — U! = U) — (jJb, with COa > Wfi. 



A. One qubit per time 

In cavity-QED, this situation is precisely the one where two- 
level atoms are sent one by one through a bimodal cavity in 
such a way that only one atom is found inside the cavity per 
time. We then let the first qubit interact with the modes. If 
the system state at i = is |OA)|OB)|ti) = |00t), where 
the modes are in the vacuum state and the qubit 1 is in the 
spin up state', we may restrict the space state to the subspace 
spanned by the basis {|10|), |01|), |00t)}. After the interac- 
tion of qubit 1 with the modes (time t), the system state will be 
given by 

m)) = ci(i)llOi) + C2(t)|0i;) + C3(i)|00t), (4) 



' We borrowed the usual notation of spin— 1/2 particles when referring to the 
states of the qubits. 



n 



ci(i) = -,^[A(1 -cos fit) +insmnt], (5) 

51^ 



C2(i) = -.SlCi(t), 

C3 (i) = -^ ( A^ + 2n'^ cos nt) , 
51^ 



(6) 
(7) 



where fl = VA^ + 217^ is the Rabi frequency. The probability 
of finding this qubit in the spin up state is 



1 



n,W = |c3WI =^(A^ + 2f7'cosm) 



(8) 



This quantity plays an important role in the generation proto- 
cols described in this section. We may use (8) to obtain the 



shortest time needed to achieve a desired value of P^ 



Ti 



A2 



t — — arccos 
from where we can see that 



A r- 

0<-<y2 



2f72 



(9) 



(10) 



1 



because the domain of the arccos function is finite. Thus, the 
qubit 1 may transfer its whole excitation energy to the field 
modes (P^^ = 0) whenever < A < \/2fi, in accordance 
with [8]. In this case, since |ci(t)p ~ |c2(t)p, the non- 
degenerate modes will end up in the Bell state 



l*+> = ^(|10) + |01)), 



(11) 



up to local unitary (l.u.) operations. It is important to remark 
that no measurements were needed here because the qubit dis- 
entangles with the modes after such t. In the cavity-QED lan- 
guage, for instance, in order to generate this state for two field 
modes, previous schemes [22, 23] generally made use of either 
two cavities (making the experimental control and costs more 
demanding) or a sudden change of the atomic frequency also in 
a bimodal setup (which can introduce errors and decrease the 
fidelity of the generation protocol). However, if one is inter- 
ested in entangling distant cavities for teleportation purposes 
[24], a two cavity setup is of course more appropriate than a 
single bimodal cavity. 

Conversely, if the first qubit ends up with P^^ > 0, it will 
be entangled with the modes. For instance, consider the case 
in which P^^ = 1/3. This may be achieved by choosing < 
A < (1 + ^/3)n and adjusting t through (9). In this case, the 
system will end up in the state (l.u.) 

|W3)-^(|io;) + |ou) + |oot», (12) 

VO 

which is a hybrid W state among these three subsystems. 
Another example is the choice P^^ = 1/2, which demands 
< A < (2 + V2)fl. Again, by adjusting t with the help of 
(9), a hybrid entangled state (l.u.) 

|M^T) = ^(|10;> + |0U)+A/2|00t)), (13) 



is generated. This is a special instance of a generalized W state, 
and it finds application in quantum teleportation [25]. 

In what concerns the generation of entangled states, the role 
of the sign of coupling constant seems to be irrelevant for the 
one qubit setup. However, when a second qubit subsequently 
interacts with the modes, the choice of different constant cou- 
pling phases will lead to substantially different dynamics. For 
completeness, we will also consider a delay time td during 
which the subsystems follow free evolution [8]. If the inter- 
action of qubit 1 with the modes takes a time ii, and the in- 
teraction of qubit 2 (initially prepared in the spin down state 
14,)) with the modes takes a time t2, the final global state can 
be written as 

|^(ii, id, t2)) - a(ii, id, t2)|10U) + /3(ii, id, i2)|01U) 

+7(ii)|00n) +<5(ii, id, i2)|00;t), (14) 

where we point out the independence of 7(ii) = c^iti) on 
subsequent times. The other coefficients are given by 

a(ii,id,i2) -ci(ii)a(t2)e-^^*'' +C2(il)6(i2)e'^*^ 
/3(ii,id,i2) =ci(ii)6(i2)e-'^*'' +C2(il)a*(i2)e^^*^ 
^(ii,id,i2) =ci(ii)ci(i2)e-*^*'' +C2(il)c2(i2)e^^*^ 

(15) 



B. N qubits at a time 

The system Hamiltonian for A^ qubits interacting with two 
bosonic modes in the rotating wave approximation reads 



N 

H = ojy^a'a + uigb'b -\ > a^ 

fe=i 

N 

+ 5]fif)(6V,-+6a+). 
fe=i 



N 

fc=i 



(aV^ 



(19) 



Please note that we are disregarding qubit-qubit interactions. 
In cavity-QED setups, for example, this corresponds to main- 
tain the atoms sufficiently apart so that direct dipole-dipole in- 
teractions may be neglected. 

We now prepare the system in the state 100)0- j^ j^)*^^, where 
each mode is in the vacuum state, qubit 1 in the spin up state, 
and any other qubit in the spin down state. We will consider 
the case where all the coupling constants have the same sign. 
Thus, similar to (2) and (3), we now move to a frame rotating 
according to 



N 



Gm = oja) a + uob'b + — N^ a^, 



(20) 



k = l 



with ci{tj) and C2{tj) obtained from (5) and (6), respectively, 
and 

a{t2) = -^[n^ + (n^ + A^) cosnt2 - iAhsinnt2],{l6) 



b{t2) = -S2-^(l - COS 17*2) • 



(17) 



We now provide a link between the sign of the coupling con- 
stants and the delay time. We will show that the choice of the 
relative sign of coupling constants may be chosen to mimic the 
delay time under special conditions. In order to see that, one 
has to realize that the probability of finding the qubit 2 in the 
spin up state P-|-2(ii, tc(,t2) ~ |(5(ii, ic(,i2)P can be written, 
for the choice td = J7r/(2A) {j = 0,1,2,.. .), as 

P^^(il,i2)-2{|ci(ii)nci(i2)P 

+ (-l)-''siS2Re[c?(ii)c2(i2)]}. (18) 

Therefore, the choice of si = — S2 together with no delay 
time (j = 0) is equivalent to si ~ S2 with j odd. The possi- 
bility of eliminating the delay time through careful choices of 
the phases of the coupling constants may be useful for dimin- 
ishing decoherence damage (reducing overall time). This is a 
quantum interference effect. 

The pursuit of generation of W states involving more qubits 
in this scenario may not be useful since the number of steps 
in this generation scheme grows with the number of qubits, 
and decoherence is expected to decrease the state fidelity if 
many steps are required. This motivates the search for schemes 
in which the generation of A^-qubit entangled states does not 
scale with N. The rest of the paper is dedicated to such 
schemes. Indeed, it turns out that the next scheme does not 
scale with N (one step generation), and the time taken to com- 
plete the protocol decreases with N. 



and the evolved state at time t will be [26, 27] 

|V'iv(i))=aAr(i)|10)|i)^^ + 5^(i)|01)|;)^^ 

N 

fc=i 



(21) 



with the coefficients 



aw(i) 



n 



[A(l — cos (1^1) + iflpf sin ^1^1] , (22) 



'AT 



n 



&w(i) = ^^[^(1 — cos Q,pft) — ifl^ smQ.pft], (23) 



2rj2 

CM'(t) = 1 - -^1^(1 - COsilNt), 



-N 



'TV 



„('=5^1) 



Jl) 



-AT '(i) = c);^(i)-i. 



(24) 



(25) 



where fijv = ^/A^~+~2NT!^. Of course, for A^ = 1 we recover 
the results of previous analysis (observe that the last equation 
above does not apply in this case). For iV > 1, it follows from 
the last two equations that 



1 Ki 

t = -^ — arccos 



P(f) - 1) 



2^(2 



202 



1 



"at 



H 



{N) 



Tfc 



2n^ 



n 



N 



2f]2 



(26) 



(N) 



is the shortest time needed to achieve a given probability P^ 

of finding qubit 1 (A: 7^ 1) in the spin up state. Also, from the 
arguments of the arccos functions above, these two probabili- 
ties are related by 



P. 



(N) 



Ti 



1 



P 



(N) 



(27) 



As we are interested in the generation of hybrid (prototype) While for a p imaginary we have 



:,{N) 



(N) 



WAr+2 (Wat) States, we must impose P|"^ = P\" ' = 1/(A^ + 
2) (P.^^^ = P.'^' = l/N) which is only fulfilled for iV = 2 

I 1 I A: 

{N — 4), see (27). Thus, for A^ > 1, this initial prepara- 
tion only allows for the generation of hybrid or prototype W 
states involving either two modes and two qubits or four qubits 
(modes factoring out). Such a limitation may also be observed 
in similar situations where only one excitation is present in the 
system [28, 29]. Here, the hybrid (prototype) W4 state may be 
generated for a detuning satisfying < A < 20 (A = 0). 

The limitation on N for generation of hybrid (prototype) W 
states may be avoided if we follow a prior step: we let an 
auxiliary spin up qubit (frequency wg) interact with the modes 
(initially in vacuum states) transferring its whole energy to the 
modes [(9) and (10) with P^^ (t) = 0]. This is the exact proto- 
col described in the last section. Conveniently parameterizing 
the resulting modes state as \ipc) ~ p|10) — soP*|01), where 
So = —1 is the coupling sign of the auxiliary qubit with mode 
Pand 



p = --{Ao/n + i 



2-A2/f^2 



(28) 



with Aq — LOA — <^Q ^ <^Q — ^B, it follows that |pp = 1/2 
whenever < Aq < \/2^, leaving the fields in \ipc) which is 
l.u. to the Bell state (11). The auxiliary qubit at the end of this 
step is in the spin down state and is disregarded thereafter. 

We now let N qubits, each of them initially in the spin down 
state, simultaneously interact with the entangled modes. We 
again consider the case where all the coupling signs Sk are pos- 
itive. The evolved state of the system after an interaction time 
t is still given by (21), except for the coefficients that now read 

aN{t)^^{{p-p*)Nn' + [pA^ + {p + p*)NQ.^] 

"at 

X cosr^Ari — ipAr^AT sinl7jvt}, (29) 

6jv(t) = ^{{p* -p)Nn^ + [p*A^ + {p + p*)Nn^] 



X cosCInI + ip* ACIn sin J7 Aft}, 

n 



(30) 






X sin JItv^]- 



(31) 



Please note that there is no differentiation between the ampli- 
tudes of probability (referred to the spin up state) for the qubit 
1 and the others. This fact will give us freedom to generate W 
states with any number of qubits. In the following we consider 
p real (Aq = \/20) or imaginary (Aq ~ 0). For a p real, it 
follows that 



np 



p 



1 "w 

t = — — arcsin 

iljv V 20 



(N) 



tk 



(32) 



is the shortest time needed to achieve a given probability 
pI ' (t) of finding qubit k in the spin up state, with a detuning 

I k 

in the interval 



O^^A 



P. 



(N) 



N. 



(33) 



1 
t = — — arccos 

OjV 



V2nA-n%Jp' 






\/20A 
with the detuning now in the interval 



(34) 



,(1-/ 



1 -NP' 



(^^)^ 



Tfc 



< 



P 



(N) 



A 

n 



< 



.(! + / 



1 -NP 



Wn 



Tfc 



P 



(N) 



(35) 
In order to generate WAr_|_2 (Wjv) states, we must impose 



{N) 



1/{N + 2) (P 



(N) 



l/N) in the above equations. 



This implies that hybrid (prototype) W states occur if the de- 
tuning is such that < A < 20 (A = 0) for p re al o r 

(\/2a/]VT2 - 2)0 < a < {^/2y/W 



2 + 2)0 (A = V27V) 
for p imaginary. Careful analysis of (32) and (34) reveals that 
the time required to generate Wn+2 or Wn states decreases 
with N . A similar result is presented in [30]. We emphasize 
that due to its importance in the multipartite entanglement the- 
ory [31], many different schemes for generation of W states 
have been previously proposed [32^0], and their experimen- 
tal realization has also been reported [41^5]. 



III. GENERATION OF W, GHZ, AND CLUSTER STATES 
WITH BIMODAL DISPERSIVE INTERACTIONS 

It is a well known fact that a single bosonic mode may medi- 
ate interactions between two previously non-interacting qubits 
[46]. Consider for example that the mode frequency is suffi- 
ciently apart from the qubit frequencies (taken to be all equal), 
the so-called dispersive limit. In the two-qubit case, the mode 
in the vacuum state then induces a dipole-dipole interaction 
between the two qubits (XY interaction), allowing the gener- 
ation of a Bell state at certain times of interaction [46]. In the 
following, we consider the dynamics far from resonance of two 
modes simultaneously interacting with N qubits, and explore 
the role played by the phases of the coupling constants. 

The simultaneous interaction of N qubits with two non- 
degenerate modes is described in the interaction picture by the 
Hamiltonian 



N 



Hi{t) = O ^(aV^e^^* + aa+e-'' 



-iAt\ 



k=l 

N 

+O^Sfc(6V^e-*^* + 6a+e''^*), (36) 

k=l 

where we have again set the qubit frequencies (cj) midway the 
modes frequency separation, so that A = uja — uj = uj — ujb, 

with W^ > UJg. 

The dispersive limit is obtained when the detuning A is suf- 
ficiently larger than the coupling constant O. Following the 
method described in [47], we obtain from (36) the effective 
Hamiltonian 



N 



N 



H,,, = -\{a^ a-h'^b)Y,al- \Y,{1- s,Sk)o+a- , (il) 



k=l 



j,k 



where A = A^/0 is the effective qubit-qubit coupling constant 
induced by the modes. 



A. W states 



B. GHZ states 



For A^ == 1, the above Hamiltonian has already been con- 
sidered for the generation of multidimensional entangled co- 
herent states of two cavity modes [48]. The first term of the 
above equation describes the Stark-shifts. These induced shifts 
may cancel each other if the modes are prepared in Fock states 
with equal number of photons (or simply in the vacuum). The 
second term (dipole-dipole interactions induced by the modes) 
may be totally inhibited if Sfe = 1 for all k. This occurs because 
the interaction between qubits j and k induced by the mode A 
cancels the interaction induced by the mode B. Clearly, this 
effect is absent when only one mode interacts with the qubits. 
We may then state a general rule that given two qubits with 
Sj and Sfe, the induced dipole-dipole interaction between them 
will be inhibited whenever Sj — s^- The contrary is also true, 
i.e. their induced dipole-dipole interaction will be doubled pro- 
vided Sj — —Sfe. 

We still consider cancellation of the Stark-shifts in (37). We 
then consider that only one of the N qubits is initially in the 
spin up state, say the first qubit, and the others are in the fun- 
damental state, so that the initial state of the qubits is a^\\.)^^ . 
If we set Si = — 1 and Sfe = 1, for fc > 1, the state of the N 
qubits after simultaneously interacting with modes will be (at 
a time t) 



N 

iV'ivW) = E4''(^Kii) 



(»N 



k = l 



with coefficients 






V^v^l 



(38) 

(39) 
(40) 



From this follows that 



arccos 



t = 



P. 



(N) 



Ti 



2^/w^^:x 



(41) 



(N) 



is the shortest time needed to achieve a given Pi , i.e. the 
probability of finding the qubit 1 in the spin up state. 
interesting to notice that now 



It is 



l-Pf)^(7V-l)pW, 



so that the qubits after a time t are left in the state 



(42) 



\^N{t)) = / 



P 



(N)„+l l\(SN 



l-E 



(N) 



ti 



i;)^ 



+ ^\ 



Ti 



N -1 



N 

E 

fc=2 



4li) 



0S>N 



(43) 

This means that a wide class of W states may be generated us- 
ing the present scheme. For example, if A^ = 3, we may gen- 
erate the usual (prototype) W3 or its variant (suited for telepor- 
tation) just by setting Pl ' = 1/3 or Pi ' — 1/2, respectively. 
As in the last section, analysis of (41) also reveals the decreas- 
ing of interaction time with the number of qubits. 



This configuration is also useful for generating more kinds 
of multipartite entangled states when we leave the one- 
excitation sector of the system Hamiltonian. For example, if 
the qubits are prepared in the state |+i) 1+2)1+3) = |+++), 
with \+j) = (1/V2)(|t^-) + ||)j), the evolved state after a time 
t will be 



mt)) 



1 



■.[\m) + ip + ^'y)\m) + (p + ^'y)\m) 



2^2 

+{^i + i,^)\m) + (/i + ii^mt) + i;u)], (44) 

where /i = cos 2\/2At and u = (l/\/2) sin2\/2At. Hence, if 
we let the qubits interact with the modes for At ~ tx l{2\fi), 
the following state will be generated 



IV') 



1 



2^2 



(int) - im) - int) - im) 



-i;n) - i;n) - lut) + lu;)), 



(45) 



which is l.u. to the GHZ state, e.g. VU\i)) = (l/\/2)(|ttt) + 
1;;;)), where U = {-^ialfl'^ ® {iolf'^ ® {i<^lV^^ and 
V — (— \/i)o'f®H2® "Ha, with "Hi representing the Hadamard 
operation on qubit i. It is known that U belongs to a set of local 
Clifford unitaries mapping the whole equivalence class associ- 
ated with a given graph state [49-5 1 ] . The GHZ state has an 
important role in the study of fundamentals of quantum me- 
chanics. Greenberger, Home and Zeilinger (GHZ) have shown 
that the predictions of local realistic theories and quantum me- 
chanics are completely different using this state [52]. Other 
theoretical proposals for generation of GHZ states have been 
previously published [53-56] as well as its experimental real- 
ization [57-61]. 



C. Cluster states 

Now consider A^ = 4, si = S2 = ^ 1 and S3 = S4 = 1. In 
this case, Hamiltonian (37) reads (modes in the vacuum) 

iJcff = -2\[{ala^ + (J^cr^) + {cri(T^ + erf cr;^) 

+ {a^(J^ + '72"cr;^) + (c^:^crr + <^2'^4)]- (46) 

If the initial state of the qubits is |t|ti)^ the evolved state at 
At = 7r/(4\/2) willbe 

w = Imu) - imt) - i;m) - i;nt)), (47) 

which is l.u. to the usual form of a linear cluster state [13], 

e.g. -af ® 12 ® fTf ® i4i^) - (i/2)(inn) + itm) + 

liitt) ^ liiii))- Multidimensional cluster states form the ba- 
sis of what is called one-way quantum computing [62]. By 
initializing the system in these highly entangled states, and by 
properly choosing local measurements on the qubits (plus clas- 
sical feed-forward), it is possible to achieve a universal model 
for quantum computing [62]. Previous schemes for generation 
of these states have been proposed [63-67], and experimental 
realization has also been reported [68-71]. 



IV. GENERATION SCHEMES IN THE PRESENCE OF 
IMPERFECTIONS 

In this section we analyze the effects of usual imperfections 
that may arise in a realistic setup. Whereas the preceding re- 
sults were put in the general language of bosonic modes and 
qubits, the sources of imperfections we are going to discuss 
are well suited for the cavity QED context, where two-level 
atoms (qubits) can interact with bosonic modes sustained by 
a cavity. We separate the analysis of imperfections in three 
subsections. In Section IV A, we investigate the effect of dis- 
sipative environments in the generation schemes of previous 
sections. In Section IV B, we estimate the error in the state 
preparation arising from the intrinsic uncontrollability of the 
time of flight of atoms. Both sources of imperfections will be 
analyzed in terms of the fidelity F = Vv/TpIW [2]> where ji/;) 
is the ideally generated (target) state and p is the result in the 
presence of imperfections. In Section IV C, we show the pos- 
sibility of violation of local realism considering our generation 
protocol of cluster states in the presence of such imperfections. 



A. Effects of dissipatlve environments 

In real experiments, one usually faces the effects of many de- 
coherence channels. For the system considered here, the most 
relevant damages to the dynamics come from photon leakage 
from cavity mirrors and fluorescent radiation from the atoms 
(spontaneous emission). The usual formalism to treat irre- 
versible decaying processes considers the coupling of the sys- 
tem to large reservoirs, generally modeled as baths of harmonic 
oscillators [72]. Here we consider zero temperature local reser- 
voirs for each subsystem. 

We start by considering dissipation in the generation 
schemes of Section II. In the rotating-wave and Born-Markov 
approximations, the density operator for the system obeys the 
following master equation [72] 

N 

p{t) = -i[Hk,p{t)] + CApit) + CBp{t) + ^ Ckp{t), 



fe=l 



where Hk is given by (3), 

^Ap{t) = '^{[a,p{t)a)] + [ap(t),a^]), 



2 

Kb , 



CBp{t) ^ ^{[b,p{t)b^ + [bp{t),b^), 

^kPit) = J(KT,p(t)a+] + KTp(t),(T,7]), 



(48) 

(49) 
(50) 
(51) 



are the Lindblad operators acting on the density operator p{t), 
i^A, i^B^ and 7 are the decay rates for cavity mode A, B and 
the atoms, respectively. 

In Figure 1, we plot the fidelity between p and j?/') as a func- 
tion of x/^ for three different target states. In these calcula- 
tions, p is also considered at an interaction time given by (9) 
with A = ^/2Vt. The parameter x is used to embody different 
configurations for the decay rates. These plots show that the 
worst situation corresponds to that of equal dissipation rates 
(cavity modes and atoms). It is noteworthy that the generation 
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FIG. 1. (Color online) Fidelity for the cavity Bell state (top), hybrid 
W3 state (middle), and hybrid Wt state (bottom). The configurations 
are x = ka = ks = 7 (black solid line), x = O.Ika = O.Iks = 7 
(blue large-dashed line), x ~ i^a = kb = Q-^'y (red dashed line), 
X — O.Ika = O.Sks = 7 or X = 0.5ka — O.Iks = 7 (green 
dot-dashed line), and x ~ i^a = 0.5kb = O.I7 or x = O.Ska = 
Kb ~ O.I7 (brown dotted line). 



protocol aiming at the Wt state (13) is more robust to dissipa- 
tlve effects when compared to the ones for W3 state (12) or the 
cavity Bell state (11). 

In order to compare the quality of the protocols for gener- 
ation of W and GHZ states under dissipation, we move now 
to the cases treated in Section III. In these cases, the master 
equation (48) is still valid with the only modification being the 
replacement of Hk for Hctf, the latter given by (37). Once the 
modes are still prepared in vacuum states, we may still disre- 
gard the Stark-shifts appearing in (37) as well as any losses 
due to cavity damping. It can now be seen from Figure 2 that 
the W states protocols are again more robust in the presence of 
losses, now compared to the protocol for generating the GHZ3 
state. Notably, the generation scheme for Wt is still the most 
robust. 

The analysis of imperfections for the cluster state is pre- 
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FIG. 2. (Color online) Fidelity for the GHZ3 state (black solid line), 
W3 state (blue large-dashed line), and Wt state (red dashed line). 



sented in a separate section later, where we discuss applica- 
tions in experiments for testing violation of local realism. 



B. Errors in the control of atomic time of flight 

In majority of the generation schemes presented in this arti- 
cle, two or more atoms are required to simultaneously interact 
with the cavity modes. In realistic setups, however, it might 
happen that the atoms will not enter the cavity exactly at the 
same time, consequently inducing an error in the target state 
generation. In the case involving only two atoms, the inclu- 
sion of such error source is a simple matter as discussed, for 
instance, in reference [46]. The approach consists of suppos- 
ing one of the atoms reaching the cavity in advance or in delay 
compared to the alleged initial interaction time. 

We need to generalize this approach in order to analyze this 
source of error in our protocols. We do this by choosing the 
time each atom enters the cavity from a normal (Gaussian) 
distribution centered at i = (the ideal time for entering the 
cavity). The uncontrollability of the time of flights translates 
then as finite standard deviations for the normal distribution. 
In Figure 3, each point in the plotted fidelities corresponds to 
an average over many realizations, each one having the same 
standard deviation measured as a fraction (percentage) of 1/A. 
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FIG. 3. (Color online) Average value of fidelity as a function of the 
fluctuation of the atomic times of flight measured as a fraction (per- 
centage) of 1/A for GHZ3 (black circles), W3 (blue triangles), and 
Wt (red squares). Here, 3, 000 repetitions have been considered to 
evaluate the averaged fidelity. 

It is a remarkable fact that the generation of Wt and W3 are 
now affected in almost the same way by moderate fluctuations 
in the interaction times. Also, the GHZ3 state protocol, which 
was previously shown to be more sensitive to energy relax- 
ation, turns out now to be more robust than the other protocols 
for errors in the time control. 



C. Cluster state nonlocallty under Imperfections 

In reference [73], Scarani, Acin, Schenck and Aspelmeyer 
introduced a nonlocality test (hereafter mnemonically referred 
to as SASA inequality) suitable for cluster states. The SASA 
operator is given by [73] 
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FIG. 4. (Color online) Average value of the SASA operator as a func- 
tion of 7/A for different levels of precision in the control of the time 
of flight of the atoms, measured as percentages of 1/A, for 0% (black 
dashed line), 5% (blue triangles), and 10% (red squares). Here, 3, 000 
repetitions have been considered to evaluate the average (B) . 



and any local realistic theory will lead to {B) < 2. For the 
cluster state 

I04> = ;^|+)|t)l+)lt) + ^l+)lt)|-)l;) 



1 



-)it) + ;,i-)i;)i+)i;>, 



(53) 



where |±) = (|t) ± |4))/\/2, it is found maximal violation, 
i.e. (B) — 4. The experimental verification of violation of the 
SASA inequality by linear cluster states l.u. to (53) has been 
demonstrated [74, 75]. 

In order to evaluate the expectation value of the SASA op- 
erator in form (52), using the cluster state 1-0) generated in our 
protocol, all one have to realize is that T\tp) — {(p^) with 



T 



-■HifTi «)l2«)Cr3 (8)^4- 



(54) 



Following the same guidelines explained in the last subsec- 
tions, we now analyze the effect of dissipation and time of 
flight fluctuations on the violation of the SASA inequality as 
shown in Figure 4. It is a remarkable fact that (B) is quite ro- 
bust against error in the control of atomic times of flight. Be- 
sides, the nonlocality of the state generated under these imper- 
fections is revealed even near the weak coupling regime where 
the dissipation constants are comparable to the coupling con- 
stants. To be more specific, (B) < 2 for decay constants only 
smaller than 7 « 0.4A. 



V. CAVITY QED SETUPS 



6 = 



I2 



I2 



+al (g) cr| (g) cr| (g) 0-4 - tJi (g) cr| (g) CTg CT^, (52) 



Since some protocols presented in this work are supported 
by the ability of willingly choosing magnitudes and phases of 



cavity 1 




FIG. 5. Sketch of the setup proposed in [8] showing cavities 1 and 2 
and the dielectric barrier. The degenerate eigenfunctions are pictori- 
ally depicted. 



atom-field coupling constants, we need now to discuss these 
assumptions on potential physical setups. 



A. Proposition I 

One possible setup consists of two weakly coupled identical 
cavities 1 and 2 and flying atoms [8]. These cavities are actu- 
ally the two halves of a perfectly conducting rectangular res- 
onator divided by an insulating barrier of suitable thickness and 
dielectric constant. The experimenter is then endowed with 
the choice about which cavity each atom will be sent through. 
The special feature of this system that allows one to willingly 
chose the coupling constant phases is that two weakly coupled 
identical cavities possess two nearly degenerate eigenfunctions 
(frequency split is approximately proportional to the amplitude 
transmissivity of the barrier): one symmetric and the other an- 
tisymmetric, with respect to the junction of the cavities, as 
shown in Figure 5. Consequently, when sent to cavity 1, an 
atom interacts with two positively-valued mode functions, and 
when sent to cavity 2 it interacts with modes having opposite 
values. This fits perfectly the needs for phase control found in 
our state generation protocols. 



B. Proposition II 

Now we consider a one-dimensional model describing the 
interaction of a two-level atom j and M field modes inside 
a cavity [76]. Such a simplified version of the general two- 
or three-dimensional model still reveals the essential features 
of the dynamics [76, 77]. Under the dipole and rotating-wave 
approximations, the interaction Hamiltonian can be cast in the 
form 



M 



Hint = Yl gi^Haicr^ + an(J+), 



(55) 



where the position-dependent coupling constant for the atom j 
and mode n is given by 



with fc„ — Wn/c = mr/L and dig denoting the dipole matrix 
element of atom j. Once we are dealing with identical atoms, 
we may suppose equal values for any atomic dipole matrix el- 
ements deg . Under appropriate choices of polarization and by 
considering the atomic frequency midway the separation of the 
natural frequencies of consecutive cavity modes, say n and 
n + I, the above Hamiltonian may be reduced to effectively 
contain two modes 

Hint = g^nKa^aJ + a„cr+) + ffi+i(4+i<^7 + an+ia+). 

(57) 
The problem consists in finding solutions for equations of the 
kind 
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U) 



UJn_ 



4n^sin(fc„rj), 



(56) 



k„ sin[fc„(fj + 1/2)] = ±y fc„+i sin[A:„+i(fj + 1/2)], 

(58) 

for any atom j and two selected modes, where we have defined 
parameters scaled by the cavity length, i.e. fc„ = Lkn = mr 
and fj — rj/L and set the origin rj = to be the center of 
the cavity, i.e. the cavity mirrors are located at — L/2 and L/2. 
From now on, we will consider the lowest order modes 1 and 
2 to illustrate the general idea. 

Let us consider first the scheme in Section II A, where one 
atom (j ~ 1) interacts with each of two modes with coupling 
constants with either equal or opposite phases. In this case, we 
must solve the simple equation 

/fci sin[fci(fi + 1/2)] = ±Y fca sin[fc2(ri + 1/2)]. (59) 

It is easy to find that ri = T^ ai'csin(— t=), with the minus 
referring to the case with equal phases and the plus the other 
way round. Of course, other positions for the atom may be 
obtained if we select modes with higher n and n + 1. 

For the schemes with two or more atoms crossing simul- 
taneously the cavity (Sec. refsec:disp), we have to find so- 
lutions of system of equations involving terms of the kind 
(58). In the case of two atoms there is a solution valid for 
equal couplings with opposite phases [f i = ^ arcsin( —K= ) 

and f2 = — :^ arcsin(— T=)]. However, in contrast with the 
first proposition, any other solutions are not physical for this 
setup . The reason is that they imply putting more than one 
atom at the same point in the cavity axis. This is a very un- 
favorable situation because, even if we tolerated some error 
in the magnitude of the coupling constant by considering both 
atoms at slightly equal positions, we would still be in trouble 
because direct dipole-dipole coupling could not be neglected. 
This would alter the Hamiltonians used in the protocols in a 
very substantial way. There might still be a solution for this 
problem in the 2D and 3D versions of the model, but we leave 
this possibility for future investigations. 



VI. CONCLUSION 

We first considered the quasi resonant interaction of two 
modes with either one qubit at time or with N qubits. In the 
one excitation sector of the system Hamiltonian, we found the 
deterministic generation of bosonic Bell states between two 
modes. W states may also be generated either among the two 
modes and the qubits or among the N qubits. However, when 



considering the simultaneous interaction of the modes with the 
qubits, we observed that the generation of W states occurs only 
for A^ = 2 or A^ = 4. To circumvent this constraint, we 
slightly modified the scheme by first preparing the modes in 
a Bell state. Thus, with the excitation initially stored in this 
entangled state, Wjv+2 and Wn states are generated for all N. 
Also, we observed that the generation of these states becomes 
faster with the increase of N. 

In the far off resonance interaction, we analyze the effective 
dynamics when two non-degenerate modes interact simultane- 
ously with N qubits. When only one excitation is present, it 
is possible to generate a wide class of W states involving N 
qubits, and again the interaction time decreases with N. We 
also observed the generation of GHZ states among three qubits 
initially prepared in coherent superpositions of spin down and 
spin up states. With two excitations in the system, we showed 
the one-step generation of cluster states among four qubits. 

We would like to remark that all the protocols for generation 
of multipartite entangled states presented here are determinis- 
tic (no measurements). This became possible due to quantum 
interferences arising from different choices of coupling con- 
stant phases. 

We also studied these generation schemes under the the pres- 
ence of important sources of imperfections such as energy re- 
laxation and fluctuations in the time of flight of the atoms 



through the cavity. The results showed good fidelities for an 
appreciable range of the parameters involved. As an applica- 
tion of our generation protocols, we investigated nonlocality 
tests for the cluster state in the presence of those imperfections. 
The results are also very robust showing violation of the S AS A 
inequality for dissipative rates very close to the weak coupling 
regime of cavity QED. The quest for quantum effects in typ- 
ically unfavorable situations (like the weak coupling regime) 
is a very active research area [78]. The violation of SASA in- 
equality under such severe imprecisions makes our protocols 
useful for this kind of investigation. 
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